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ABSTRACT 

In this paper we develop the extended jet model of Potter & Cotter to model the simultane- 
ous multi-wavelength spectra of six Compton-dominant blazars. We include an accelerating 
parabolic base transitioning to a slowly decelerating conical jet with a geometry set by obser- 
vations of M87 and consistent with simulations and theory. We investigate several jet models 
and find that the optically thick to thin synchrotron break in the radio spectrum requires the 
jet to first come into equipartition at large distances along the jet, consistent with the observed 
transition from parabolic to conical at 10 5 i? s in the jet of M87. We confirm this result analyti- 
cally and calculate the expected frequency core-shift relations for the models under consider- 
ation. We find that a parabolic jet transitioning to a ballistic conical jet at 10 5 i? s , which starts 
in equipartition and becomes more particle dominated at larger distances, fits the multiwave- 
length data of the six blazars well, whilst an adiabatic equipartition conical section requires 
very large bulk Lorentz factors to reproduce the Compton-dominance of the blazars. 

We find that all these blazars require high power (> 10 39 W), high bulk Lorentz factor 
(> 20) jets observed close to the line of sight (< 2°) as we expect from the blazar sequence 
and consistent with the results from Paper II. The inverse-Compton emission in our fits is due 
to inverse-Compton scattering of high-redshift CMB photons at large distances along the jet 
due to the high bulk Lorentz factors of the jets. We postulate a new interpretation of the blazar 
sequence based on the radius of the transition region of the jet (where the jet is brightest in 
synchrotron emission) scaling linearly with black hole mass. 

Key words: Galaxies: jets, galaxies: active, radiation mechanisms: non-thermal, radio contin- 
uum: galaxies, gamma-rays: theory, black hole physics. 



1 INTRODUCTION 

Blazars are the most luminous active galactic nuclei (AGN). Super- 
luminal motion is observed in their jets indicating that the jets are 
relativistic and observed at small angles to the line of sight. The 
emission from the relativistic jets is strongly Doppler-boosted so 
that the non-thermal radiation overwhelms the observed emission 
from the host galaxy. High energy electrons in the plasma jets emit 
synchrotron and inverse-Compton radiation resulting in the charac- 
teristic double-peaked spectrum, with synchrotron radiation from 
radio to UV/x-rays and inverse-Compton emission from x-rays to 
very high energy 7-rays. Blazars are observed to have self-absorbed 
synchrotron radio spectra which are nearly flat in flux. 

In the first two papers in this series l |Potter & Cmter| ( (20T2t 
and |Potter & Cotter] < |2013| >, hereafter Paper I and Paper II respec- 
tively) we developed an extended jet model based on recent VLBI 
observations of M87, numerical simulations and theory. In Paper I 
we introduced a simple ballistic conical jet model which we showed 
could successfully reproduce the quiescent spectrum of BL Lacer- 
tae across all wavelengths using a relatively low power, low bulk 
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Lorentz factor jet. We included a thorough treatment of emission 
processes by integrating the line of sight synchrotron optical depth 
through the jet and numerically integrating the exact expression 
for synchrotron self-Compton (SSC) emission. In Paper II we de- 
veloped this model to include an accelerating parabolic base tran- 
sitioning to a slowly decelerating conical jet motivated by radio 
observations of M87 (Asada & Nakamura (201 l| l). We included 
inverse-Compton emission from synchrotron, CMB, accretion disc, 
starlight, broad line region (BLR), dusty torus and narrow line re- 
gion (NLR) scattered seed photons by Lorentz transforming the 
photon distribution into the plasma rest frame and using the full 
Klein-Nishina cross-section. 

We used this model to successfully fit to the quiescent spec- 
trum of the most Compton-dominant object in the FERMI sample 
PKS0227-369 ( |Abdo et al. | pOTO)), using a jet geom etry fixed by 
the radio observations of M87 (Asada & Nakamura (20 Hi) scaled 
linearly with black hole mass. We found that the fit to PKS0227 
required a high power, high bulk Lorentz factor jet observed close 
to the line of sight, consistent with our expectations from the blazar 
sequence (Fossati et al.| (1998) ) and unification ( Urry & Pa dovani] 
<jT^95)). We found that the spectrum of PKS0227 was well fit- 
ted by a jet which comes into equipartition outside of the BLR; 
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within the dusty torus or at a larger distance along the jet. We found 
that the inferred black hole mass of our fit was close to a previ- 
ous estimate using a transition region (where the jet comes into 
equipartition and the jet transitions from parabolic to conical) at 
a distance of 10° R s , as observed in the jet of M87. For this fit 
the transition region occured at a distance 34pc from the central 
black hole where the Compton-dominance of the blazar was due to 
inverse-Compton scattering of high-redshift CMB photons which 
were strongly Doppler-boosted into the plasma rest frame due to 
the jet's high bulk Lorentz factor. We postulated that one of the rea- 
sons Compton-dominant blazars tend to be observed at large red- 
shifts is due to the temperature dependence of the CMB on redshift 

OcMB OC (l + z)). 

In this paper we set out to investigate whether this model is ca- 
pable of fitting to the spectra of a sample of Compton-dominant 
blazars and to investigate what physical properties these blazars 
have in common. We wish to find out whether the physical param- 
eters we found for PKS0227 in Paper II hold for other Compton- 
dominant blazars, in particular, if we can distinguish between mod- 
els in which the transition region occurs within the dusty torus or 
further along the jet using the optically thick to thin synchrotron 
break. This is the first time a realistic extended jet model has been 
used to fit to the spectra of a sample of blazars, so it is interesting 
to see whether our model, which is based on the observed geometry 
of the jet in M87, is capable of fitting to a sample of real blazars. 
For this sample we choose the six most Compton-dominant blazars 
(after PKS0227) from Abdo et al.| ( |2010"l > in order to test our model 
and predictions from Paper II. 

Compton-dominant blazars are thought to represent the most 
intrinsically powerful jets in the blazar sequence (Fossat i et al.| 
( |1998[ >). In this scenario we use the radio luminosity of blazars 
as a proxy for their kinetic luminosity and take an average of the 
multi-wavelength spectra in bins of radio power. A general pro- 
gression towards higher Compton-dominance and decreasing peak 
frequency of synchrotron and inverse-Compton emission for more 
powerful objects was found by Fossati et al. ( 1998}. It is currently 
believed that the Compton-dominance of powerful blazars origi- 
nates from inverse-Compton scattering of external photons ( Meyer 
|et al.|p012| >). It is thought that powerful objects with more lumi- 
nous accretion discs and broad line regions (BLRs) create a larger 
seed photon distribution for inverse-Compton scattering than in BL 
Lac type objects where the dominant photon source is synchrotron 
seed photons. It is believed that radiative cooling of accelerating 
electrons via inverse-Compton scattering of BLR photons leads to 
the lower inverse-Compton peak frequencies in Compton-dominant 
objects jMaraschi et al.| ( |2008[ >). 

In Paper II we considered scattering of synchrotron, accre- 
tion disc, BLR, dusty torus, NLR, starlight and CMB photons and 
found that PKS0227 was best fitted by inverse-Compton scattering 
of CMB or dusty torus seed photons. We found that the low inverse- 
Compton peak frequency was due to the lower peak energy of CMB 
and dusty torus photons relative to synchrotron seed photons. We 
found that the Compton-dominance was due to the increase in exter- 
nal photon density due to Doppler-boosting in the high bulk Lorentz 
factor jet. The low synchrotron peak frequency was due to the base 
of the conical section of the jet having a lower magnetic field due 
to the high black hole mass and the assumed linear scaling of the 
M87 jet geometry with black hole mass. This was the first time that 
such a scenario has been suggested for the Compton-dominance of 
blazars, however, this was also the first time that a realistic model of 
an extended jet had been used to model the spectrum of a Compton- 
dominant blazar. 



In this paper we will first briefly introduce our jet model. 
We compare several adiabatic and ballistic jet models to the radio 
emission of the six most Compton-dominant blazars in Abd o et aT] 
p0T0l > (after PKS0227). We use the slope of the radio emission 
and the frequency of the optically thick to thin synchrotron break to 
distinguish between a transition region within the dusty torus or at 
futher distances along the jet. We then show the results of fitting the 
surviving models to the multi-wavelength spectra of the six blazars. 
We compare the physical parameters inferred from the model fits to 
the results of Paper II in order to determine whether our model is 
capable of reproducing the quiescent spectra of a sample of blazars. 
We comment on our results in light of current ideas on the blazar 
population and unification. 



2 JET MODEL 

To fit to the spectra of the Compton-dominant blazars we use a 
model with an accelerating parabolic base transitioning to a slowly 
decelerating conical jet with a geometry set by the recent radio ob- 
servations of the jet of M87 ( Asada & Nakamura ( 201 It). We will 
briefly describe the principles of this model in the following section 
(for full details see Papers I and II). 

Our model is motivated by and consistent with observations, 
simulations and theory. The jet is assumed to start magnetically 
dominated at the base. The magnetic energy is converted into bulk 
kinetic energy in the plasma via a magnetic pressure gradient. 
The plasma is accelerated until it approaches equipartition between 
magnetic and particle energies where the jet approaches a termi- 
nal bulk Lorentz factor and transitions from parabolic to conical 
( |McKin ney ( 2006 1). In the conical section the jet slowly decelerates 
due to interactions with its environment (Laing & Bridle (2002), 
|Hardcastle et al.| ( |2005| > and |De Young] (2010[ l) and this converts 
bulk kinetic energy into accelerating electrons via shocks. Deceler- 
ation is observed along jets and this is also consistent with evidence 
for in situ acceleration of electrons from observations of optical 
synchrotron emission (Meisenhei mer et al.| ( |1997) and |Jester et al.| 
( |2001[ >), We assume that the plasma is isotropic and homogeneous 
in the rest frame of each section of plasma so can be described by a 
relativistic perfect fluid. We conserve energy-momentum and elec- 
tron number via the conservation equations 

V ll T^(x)=0, V^(x) = 0. (1) 

We calculate the synchrotron emission and opacity in each 
section by Lorentz transforming into the instantaneous plasma rest 
frame. We integrate the synchrotron optical depth through the jet to 
each section to calculate the observed synchrotron emission. In the 
model we treat inverse-Compton scattering of synchrotron, CMB, 
starlight, accretion disc, BLR, dusty torus and NLR seed photons 
by Lorentz transforming into the plasma rest frame and using the 
full Klein-Nishina cross-section. We evolve the electron population 
along the jet taking into account radiative losses from synchrotron 
and inverse-Compton emission and adiabatic losses due to expan- 
sion of the jet. We set the geometry of our model using the obser- 
vations of the geometry of M87 with lengths scaled linearly with 
black hole mass. 

2.1 Model parameters 

We consider a jet with total lab frame power Wj, jet length L and 
conical half-opening angle ^'opening observed at an angle ^observe 
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Figure 1. This figure shows a schematic diagram of our jet model. 

to the jet axis. We assume that electrons accelerated due to shocks 
(or via magnetic reconnection) along the jet have an initial distribu- 
tion, N e (E e ) oc E~ OL e~ E ' 11 / - Bmax , with a minimum energy E m i n . 
The base of the jet has an initial bulk Lorentz factor 70, which ac- 
celerates to 7 m ax at the end of the parabolic section. The jet slowly 
decelerates from a bulk Lorentz factor of 7 max to 7 m in by the end 
of the conical section. The central black hole has a mass M. A 
schematic diagram of the jet is shown in Figure[T] 

2.2 Electron energy losses due to emission 

In this paper we have modified our treatment of electron energy 
losses from Paper II to include the decrease of individual electron 
energies with radiation losses. In Paper II we included radiative en- 
ergy losses to the electron population due to both synchrotron and 
inverse-Compton emission, however, we neglected the decrease of 
electron energy due to radiative losses and used the approximation 
that electrons radiate the majority of their emission close to their 
initial electron energy. This approximation resulted in an increase 
in speed of our jet code and did not produce a significant change in 
the emitted spectrum of the jet, however, the value of _E max is un- 
derestimated for a fit to a spectrum. To include the transfer of elec- 
trons down through energy bins as they radiate we modify Equation 
59 from Paper II to become 

AT ,c j /n AT ,c \ Ptot(x, dx',E e ) x lm 

N e (E e ,x + TWkdz ) = N e (E e ,x) - cEe(1 _ 1/dEe) 

P to t(x,dx' ,E e dE e ) x lm 
cE e dE e (l - iJ&Ee) ' 

where N e (E e , x + 7buikdx') is the number of electrons of en- 
ergy E e contained in a slab of width lm (in the x-direction) in the 
plasma rest frame, Ptot{x, dx' , E e ) is the total power radiated by 
electrons of energy E e through synchrotron and inverse-Compton 
emission in the jet section located at a lab frame distance x with a 
width da;' in the plasma rest frame and dE e is the fractional dif- 
ference in energies between adjacent electron energy bins in our 
code. We have also corrected the typo in Equations 10 and 59 of 
Paper II and changed Is to lm. The first term in Equation [^corre- 
sponds to the number of electrons lost from the energy bin through 
radiative energy losses when crossing the section of width dx' in 
the plasma rest frame. The second term corresponds to the increase 
in the number of electrons in the energy bin due to higher energy 



electrons cascading down in energy following radiative losses. If 
choosing to take into account the decrease in electron energy as the 
electron radiates, Equation 10 in Paper II should also be modified 
to include the electron gain and loss terms above. 



3 RESULTS 

In this section we shall try to constrain the properties of the 
jet by fitting our model by eye to the simultaneous multi- 
wavelength observations of the six most Compton-dominant blazars 
(after PKS0227) J0349.8-2102, J0457.1-2325, J0531.0+1331, 
J0730.4-1142, J1504.4+1030 and J1522.2+3143 from [Abdol 
[etaLlpOTO) . First we compare the radio observations with the cal- 
culated radio emission from several different jet models. We then 
use the optically thick to thin synchrotron break to attempt to con- 
strain the location of the transition region of the jet. We calculate an 
analytic formula for the observed synchrotron break frequency of a 
general jet plasma described by a relativistic perfect fluid. We use 
this to confirm the results of our model and to infer the properties 
of parabolic and conical jets by comparing our calculations with 
observations of the frequency dependent core-shift in jets. Finally, 
we fit the multiwavelength spectra using the jet models compatible 
with the radio observations and optically thick to thin synchrotron 
break. Throughout this paper we assume a standard ACDM cos- 
mology with Ho — 71kms _1 Mpc _1 and Q.a = 0.73. 

3.1 The different models 

In this Paper we wish to constrain the physical properties of the jet 
by fitting different models to the observations of the six blazars. 
In Paper II we found that a ballistic equipartition jet best fitted 
the observations of PKS0227, whilst an adiabatic equipartition jet 
struggled to reproduce its Compton-dominance with plausible bulk 
Lorentz factors. Here we will investigate models in which the coni- 
cal section of the jet is adiabatic or ballistic with a transition region 
located outside of the BLR. In these models we will also consider 
cases where the plasma deviates from equipartition along the coni- 
cal section, however, we still assume that the transition region is in 
equipartition. 



4 RADIO EMISSION 

In Paper II we found that the spectrum of PKS0227 was best fit- 
ted by models in which the conical section of the jet was ballis- 
tic. In Figure |2ji we show the radio emission from three models 
compared to the radio data of J0457. We find that the radio obser- 
vations of these six Compton-dominant blazars have radio slopes 
which are in general steeper than the radio emission produced by 
the ballistic equipartition jet model from Paper II, as shown in Fig- 
ure]^. We find that the ballistic equipartition jet tends to overpro- 
duce low frequency radio emission compared to the data. This sug- 
gests that if the jet is ballistic and conical it becomes particle domi- 
nated outwards from the transition region since this reduces the ra- 
dio synchrotron emission at larger distances, whilst maintaining the 
Compton-dominance from scattering external photons. The radio 
emission from both the adiabatic jet in equipartition and the ballistic 
jet which becomes particle dominated fit the radio data well, whilst 
the ballistic equipartition model does not. Radio observations are 
usually neglected in most previous investigations modelling blazar 
spectra. These results, however, show that it is useful to consider 
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(a) Radio emission of J0457 (b) Optically thick to thin synchrotron break in J053 1 

Figure 2. Figure a shows the results of fitting different jet models to the radio emission of J0457. We see that the ballistic equipartition jet does not match the 
radio observations well at low frequencies. The adiabatic equipartition model fits the optically thick radio slope well as does a ballistic model which becomes 
particle dominated beyond the transition region (Ue/Ug oc x~ x from xt < x > IOOiet)- In Figure b we show fits to the radio emission of J0531 using 
adiabatic and ballistic jet models with a transition region located at x^ within the dusty torus and at larger distances along the jet (as implied by a transition 
region at W 5 R$ m M87). The change in gradient of the flux from optically thick to optically thin at 3 X 10 — 5 eV is well matched by transition regions at 
large distances along the jet ~ lOOpc and excludes a transition region within the dusty torus at a few pc since the jet remains optically thick up to much higher 
frequencies than the observed break. 



the radio emission when constraining the physical properties of 
jets. In Paper II we struggled to reproduce the observed Compton- 
dominance of PKS0227 using an adiabatic conical jet in equiparti- 
tion without using very high bulk Lorentz factors. This suggests that 
the conical section of the jet may be ballistic starting in equiparti- 
tion at the transition region and becoming more particle dominated 
further down the jet as magnetic energy continues to dissipate and 
accelerate electrons. 

We find that in order to reproduce the observed radio slope of 
these blazars the equipartition fraction of the ballistic jet should 
decrease roughly as A CCltl i — Ub/U e = 2xt/(xt + x) for 
xt > x > IOOxt (where xt is the distance of the transition re- 
gion) as shown in Figure [2^. The radio data does not constrain the 
emission from the jet at distances beyond ~10kpc and so we find 
that the ballistic jet becomes particle dominated by a factor of 50 at 
a distance of IOOxt but is not well constrained after this. The bal- 
listic jet could stay particle dominated, tend back towards equipar- 
tition or become even more particle dominated after this distance. 



5 OPTICALLY THICK TO THIN SYNCHROTRON 
BREAK 

In Paper II we found that models in which the transition region ei- 
ther occured within the dusty torus or at further distances along the 
jet were able to fit well to the spectrum of PKS0227. The main dif- 
ference between the two models was the location of the optically 
thick to thin synchrotron break due to the different radii and mag- 
netic field strengths of the transition regions in the two fits. Figure 
[2|> shows the fits to the spectrum of J0531 in which the jet has a 
transition region located within the dusty torus and further along 
the jet corresponding to a distance 10 J -R 3 as in the jet of M87. We 
see that in order to reproduce the relatively low frequency of the 
optically thick to thin synchrotron break at w 3 x 10 eV, requires 
a transition region outside the dusty torus at a distance of several 
lOOpc. If the jet comes into equipartition within the dusty torus then 



the jet overproduces the radio emission due to the optically thick to 
thin synchrotron break occuring at a higher frequency than indi- 
cated by the observations. This result holds for both the ballistic jet 
which becomes particle dominated and the adiabatic equipartition 
jet. This is interesting since for these large distances the inferred 
black hole mass of the jet agrees with the black hole mass found in- 
dependently from fitting the accretion disc spectrum to J0531 only 
if the transition region occurs at 10 R s , as we have assumed based 
on observations of M87. It is worth noting that we find it difficult 
to reproduce both the radio and gamma-ray data using a jet model 
with a transition region within the dusty torus. 

The frequency at which the synchrotron emission becomes op- 
tically thin can be calculated analytically in the case of a jet with 
a total jet power Wj, bulk Lorentz factor 7buik, radius R, and 
equipartition fraction A cqu i — Ub/UeAu this case the synchrotron 
opacity can be calculated using the formula from Paper I, assuming 
the plasma can be approximately described as a relativistic perfect 
fluid in the plasma rest frame (see Paper II) and the synchrotron 
emission is approximated to occur at the synchrotron critical fre- 
quency (see Paper I). 



3uC 



2e l/2j,5/2^ 



Et = tv, 



OtAB £ 2 



3eB 



(3) 



(4) 



We can parameterise A, which is the normalisation constant 
for the electron distribution in a slab of width lm (N e = AE~ a ), in 
terms of the power law index of the electron distribution a and the 
maximum and minimum electron energies E m i n and E max . Simi- 
larly we can express the magnetic field strength in terms of other 
known quantities. 
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3Wjfi A a 



(5) 



black hole mass of our fit, 2.5 x 10 M©, and reinforces the con- 
clusion that in the case of powerful, Compton-dominant blazars the 
majority of the synchrotron emission comes from relatively large 
distances from the central black hole. 



A = 



47{U C ( 1 + A - 



E n 



(2-q) 
-El 



(6) 



In the case of a = 2 the factor (2 - a)/(C° - £w) 
should be replaced by 1/ ln(£' max / i? m i n ) in the above and follow- 
ing Equations. We can now substitute these values into Equations[3] 
and|4]to find an expression for the synchrotron opacity. 
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This expression is quite general and is valid for any homoge- 
neous plasma which is described as a relativistic perfect fluid in 
the plasma rest frame. In our present discussion we are interested in 
calculating the jet radius at which a given synchrotron frequency be- 
comes optically thick, for a plasma in equipartition. We assume that 
the plasma becomes optically thick to a frequency when fc„dx = 1, 
where both v and dx are measured in the plasma rest frame. The 
lifetime of electrons emitting synchrotron radiation at radio fre- 
quencies are long, so the length scale dx will be determined such 
that the jet radius and magnetic field strength change significantly, 
reducing the synchrotron opacity. We take this length scale to be 
the distance over which the jet radius doubles. For a conical jet 
this requirement corresponds roughly to dx — R because in the 
lab frame the conical half opening angle is expected to be approx- 
imately 9 = 1/Tbulk an d the lab frame length is related to that in 
the plasma rest frame by a length contraction dx' = 7buikdx. Cal- 
culating k v R — 1 and rearranging we find 



R = 

Uirm 3 e c 4 
V 3e 



(7T 
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87T/i mic 2 (l + Aequi) ' E, 



qui 



2-a 
max 
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6 + Q 
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3 

6 + a 

Tbulk v 



(8) 



We see that for a conical jet in which the electron distribu- 
tion, equipartition fraction, bulk Lorentz factor and jet power re- 
main constant, the radius of the jet at which the synchrotron emis- 
sion at frequency v becomes optically thick is simply inversely pro- 
portional to the frequency v, R oc v . 

Let us now use this formula to estimate the radius of the tran- 
sition region from the spectrum of J0531. We shall use the physical 
parameters from Table [T] derived from the ballistic particle dom- 
inated jet fit to 10531 shown in Figures |2jj and [3] The observed 
frequency at which the synchrotron break occurs is approximately 
3 x 10~ 5 eV = 7.25Ghz and this corresponds to a plasma rest 
frame frequency p ^ s w 27b u iki / /(l + z). In this case we find 
R fa 7.7pc. In our model using a scaled M87 jet geometry the ra- 
dius of the transition region is 2000i? s at a distance 10 5 7? s from 
the central black hole. This places the transition region at a distance 
380pc from the black hole corresponding to a black hole mass of 
~ 3.8 x 1O 1O M0. This estimate roughly agrees with the inferred 



5.1 Frequency core-shift 

We can use the general result in Equation[7]to compare to observa- 
tions of the frequency dependent core-shift in AGN jets. Previous 
calculations of the core-shift have assumed a constant bulk Lorentz 
factor (Konigl ( 1 981} and |Lobanov| (|T998 1). Observations seem to 
show that generally awe oc v~ bB , with a mean value k fa 1 and 
a tail towards lower values of k (Sokolo vsky et al7| < |201 l[ l). We 
shall use our formulae to calculate the expected frequency core- 
shift in the case of the three conical jet models in Figure |2ji: the 
adiabatic equipartition jet, the ballistic equipartition jet and the bal- 
listic particle dominated jet, and in the case of a parabolic accelerat- 
ing jet. Core-shift measurements generally probe frequencies from 
l-40Ghz which correspond to large distances (fa 2000pc for the fit 
to J053 1) down to several tens of Swarzchild radii close to the base 
of the jet (43Ghz corresponds to 14 - 23R S in M87 |Hada et al] 
\201 1| >). With this large range of distances it is interesting to calcu- 
late the expected core-shift relation for both the outer conical and 
inner parabolic sections of the jet. Assuming that the electron dis- 
tribution remains fairly constant along the region of interest we use 
the formula. 



fc„dx oc 



a: 



W- 



dx 



(1 + Aqui)- 



(9) 



R- 



"7buik^ b ; 



Let us assume that in the lab frame R oc x n and the bulk 
Lorentz factor 7buik oc x m . In this case the rest frame distance 
dx over which the jet radius changes by a fixed fraction is pro- 
portional to dx oc x/7buik oc -R 1//?l /7buik- To be compatible 
with the core-shift measurements (x oc v~ bs ), in general, requires 
that the equipartition fraction of the jet also depends on distance 
j4 cqu i oc x w . Let us also allow the remaining jet power in a section 
to change along the jet, Wj oc x u , due to adiabatic losses for ex- 
ample. We can use Equation[9]to calculate the expected core-shift 
properties of a blazar in the simplifying cases of A cqu i 2> 1, con- 
stant A cqu i and A c „ ni < 1. 



2(4 + Q ) A .»! 

(6 + o0(2n-u) + 4(™ + 2m- 1) ' cq ' 

(10) 

, A C€nl i = constant, 



k = 



2(4 + , 



(6 + a)(2ra-u) + 4(2m - 1) 



2(4 + a) 



(6 + a)(2n - u) - (2 + a)w + 8m - 



(11) 



(12) 

To calculate the core-shift relation in the case of an acceler- 
ating parabolic jet (n — 0.5, m = 0.5) we assume that the jet is 
magnetically dominated close to the base of the jet so A cqu i 3> 1, 
most of the jet power is not radiated away until larger distances, 
u — 0, and a w 2. In order for the parabolic jet to produce the 
core-shift relation k — 1 (compatible with the innermost part of 
the jet in M87 Hada et al.| |20TT|), we require that w = 1 which 
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means the equipartition fraction increases along the parabolic sec- 
tion yloqui oc x. The jet plasma actually becomes more magneti- 
cally dominated further along the parabolic section and so contains 
a smaller proportion of its energy in non-thermal electrons. Simu- 
lations find that jets start with a parabolic shape R oc x 1 ^ 2 and an 
accelerating bulk Lorentz factor 7buik oc x 1 ^ 2 (for example 



McK- 



|inney| ( [2006| >), so this result is interesting as this is the first time that 
the frequency core-shift has been calculated for a parabolic accel- 
erating jet. In our formula we have assumed that the jet is oriented 
close to our line of sight so the Doppler factor, fcoppior oc 7b u ik, 
however, this is only valid for blazars and needs to be modified for 
jets which are oriented at larger angles to our line of sight, such as 
M87. 

Let us now consider a conical ballistic jet (n = 1, m = 0) 
which starts particle dominated and becomes more particle domi- 
nated (w = —1) further along the jet. Again we assume that the 
majority of the jet power is not radiated away until further along 
the jet u = 0. In this case we assume that A cqu i <C 1. So in the 
case of the ballistic particle dominated model with a ~ 2, we cal- 



culate a core-shift relation x oc ^ obs °. In the study by Sokolovsky 



|et al.| ( [20T7) they found that the mean core-shift relation was k = 1, 
however, this distribution is relatively wide with a median value of 
k = 1.2 and a range 0.67 > k > 2 (with one outlying point at 
k = 0.29), so it does not exclude this model. Let us now consider 
the conical ballistic jet in equipartition (n = 1, w = 0) which 
suffers no electron energy losses (u = 0) and has a constant bulk 
Lorentz factor (m = 0). The early calculation by Konigll f"l981| > 
found a relation x oc v~ x independent of a for this type of jet 
model, which is consistent with observations. From Equation[8]we 
see that using our assumptions and R oc x we recover the well 
known x oc v~ bg relation. 

Finally, we consider the conical adiabatic equipartition jet with 
constant bulk Lorentz factor (n — 1, m = 0, w = 0). In this 
case the jet power decreases due to adiabatic losses on the elec- 
trons (since the jet is constantly in equipartition) so approximately 
Wj oc R- 2/3 mdu = —2/3. For the case a ~ 2 we find k = 0.69, 
which is just within the range of values of k found by Soko lovsky] 
|et al.| |20TTJ. The model which best fits the core-shift measure- 
ments is the simple ballistic conical jet in equipartition. However, 
this model does not fit the radio slopes of these Compton-dominant 
blazars well, as shown in Figure [2^. Both the adiabatic equiparti- 
tion jet and ballistic particle dominated jets fit the radio slope well 
but differ from the mean observed core-shift relation x oc u~^ s . The 
core-shift relations we have calculated for these models are both 
lower than the mean k — 1 value with both models just within 
the observed distribution of k from a sample of 20 objects from 
|Sokolo vsky et aL| ( [201 1^ . These results show that the frequency 
core-shift measurement is useful to distinguish between different 
jet models and geometries. 



6 FITTING THE MULTIWAVELENGTH SPECTRA 

In the previous two subsections we have shown that the low- 
frequency observations indicate that the jet first comes into equipar- 
tition at large distances along the jet outside of the BLR and dusty 
torus. The radio observations indicate that if the conical section of 
the jet is ballistic then it becomes particle dominated further along 
the jet and if the jet is adiabatic then the jet remains close to equipar- 
tition. We shall now attempt to fit the six multiwavelength spectra 
using models in which the jet transitions from parabolic to coni- 
cal at 10 5 i? s with either a ballistic conical section which starts in 



equipartition at the transition region and becomes particle domi- 
nated further along the jet or an adiabatic conical jet in equiparti- 
tion. 



6.1 The ballistic model 

Figure [3] shows the results of fitting the ballistic model to the six 
Compton-dominant FSRQs. The model starts in equipartition and 
becomes particle dominated at larger distances along the jet. The 
equipartition fraction changes as Ub/Us = 2xt/(xt + x) for 
x T > x > 100xt and Ub/U £ = 1/50 for x > 100:et, where 
xt is the distance of the transition region along the jet in the lab 
frame. The model fits very well to the observed spectra across all 
wavelengths. We show the physical parameters of the fits to the 
spectra in Table [l] We find that for all six blazars their inverse- 
Compton emission and Compton-dominance is well fitted by scat- 
tering of Doppler-boosted CMB photons. The different components 
contributing to the inverse-Compton emission are shown in Figure 
[4^. This is a very interesting result since it is widely believed that 
Compton-dominance is due to inverse-Compton scattering of BLR 
scattered photons (Maraschi et al. (20081). This is the first time that 
a realistic extended jet model has been used to fit to the spectrum 
of a sample of blazars and these new results show that it is impor- 
tant to consider the jet as a whole, including radio observations, to 
determine the physics governing these objects. 

We find consistently for all these blazars that fits to their spec- 
tra require high power, high bulk Lorentz factor jets observed close 
to the line of sight. This is consistent with our expectations from the 
blazar sequence (in radio power) and from AGN unification since 
we expect Compton-dominant blazars to be high power, high bulk 
Lorentz factor FRII-type jets observed close to the line of sight. We 
also find that for these objects the black hole mass (which com- 
pletely determines the geometry of the parabolic region and the 
transition region) that is inferred from the fits to the spectra are all 
large > 1O 9 M0. This is consistent with Compton-dominant blazars 
being large elliptical galaxies with powerful FRII type jets. 

Interestingly, the black hole mass inferred from the distance of 
the transition region agrees with that inferred independently from 
fitting the accretion disc spectrum to J0531 only if the transition 
region occurs at 10 5 i? s as in the jet of M87. This supports the idea 
that the transition region occurs at 10° R s in jets and scales linearly 
with black hole mass. 



6.2 The adiabatic model 

Figure [5] shows the results of fitting the adiabatic jet model in 
equipartition to the observations of the six Compton-dominant 
blazars. The model fits the data well across all wavelengths, how- 
ever, to fit the Compton-dominant spectra requires very high bulk 
Lorentz factors as shown in Table [2] The high bulk Lorentz fac- 
tors are required in order that the rate of inverse-Compton emission 
from scattering CMB and NLR photons is not much smaller than 
the rate of adiabatic losses. In the case of J1522 we find that the 
adiabatic equipartition model struggles to reproduce the Compton- 
dominance of the spectrum and requires an extreme value for the 
bulk Lorentz factor at the transition region 7b u ik = 100, this is 
much higher than the expected value of the bulk Lorentz factor in 
blazars inferred from superluminal motion Jorstad et al. (200l](. For 
this reason and the low value of k calculated in Section 5.1 for the 
frequency core-shift we disfavour the adiabatic equipartition jet. 
We find that the fits all require high power, high bulk Lorentz 
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Figure 3. This figure shows the results of fitting the ballistic jet model to the SEDs of J0349, J0457, J0531, J0730, J1504 and J1522. The model fits the 
observations very well for all six blazars across all wavelengths. We find that the inverse-Compton emission of all the blazars is well described using CMB 
seed photons with a small contribution to the highest energy emission by scattering of NLR photons. 



factor jets observed close to the line of sight as we expect for 
Compton-dominant blazars and in agreement with the ballistic fits. 
We find that all the fits have transition regions at large distances out- 
side the BLR and dusty torus with large inferred black hole masses 
> 1O 9 M . The inverse-Compton emission is due primarily to the 



scattering of NLR and CMB photons at large distances along the jet. 
We show the different components of the inverse-Compton emis- 
sion for the fit to J 1504 in Figure |4j). As for the ballistic model, 
the black hole mass inferred from the distance of the transition re- 
gion agrees with that inferred independently from the accretion disc 
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Figure 4. This figure shows the different components of the inverse-Compton emission for the fit to J1504 for the two models. Figure a shows that the inverse- 
Compton emission in the ballistic model is primarily due to scattering of CMB seed photons with a small contribution from NLR photons required to fit to 
the highest energy data point. Figure b shows that in the adiabatic model the inverse-Compton emission is primarily due to scattering NLR photons with a 
contribution from scattering CMB photons at x-ray energies. For these Compton-dominant blazars we find that the SSC emission is sub-dominant due to the 
comparitively low magnetic field strength at the transition region of the jet. 



Parameter 


J0349.8 -2102 


J0457.1-2325 


J053 1.0+1331 


J0730.4-1142 


J1504.4+1030 


J1522.2+3143 


Wi 


5.0 x 10 39 W 


1.7 x 10 39 W 


5.0 x 10 39 W 


4.5 x 10 39 W 


5.5 x 10 39 W 


1.5 x 10 39 W 


L 


5 x 10 20 m 


1 x 10 21 m 


2 x 10 21 m 


1 x 10 21 m 


1 x 10 21 m 


1 x 10 21 m 


R ■ 


5.11 MeV 


5.11 MeV 


15.3 MeV 


15.3 MeV 


15.3 MeV 


15.3 MeV 




8.10 GeV 


8.48 GeV 


4.06 GeV 


4.55 GeV 


17.3GeV 


5.48 GeV 


a 


1.5 


1.85 


1.65 


1.7 


1.85 


1.65 


9' 

opening 


3° 


3° 


3° 


3° 


3° 


3° 


^observe 


1.0° 


1.0° 


1.0° 


1.0° 


1.0° 


1.0° 


70 


4 


4 


4 


4 


4 


4 


""/max 


30 


35 


25 


28 


45 


40 


Tmin 


15 


20 


15 


15 


20 


20 


M 


8.0 x 10 9 M Q 


1.4 x 10 10 M Q 


2.5 x 10 10 M© 


2.0 x 1O 1O M 


2.5 x 10 10 M© 


1.3 x 10 10 M© 




8.3 x 10 38 W 


4.2 x 10 38 W 


3.3 x 10 39 W 


2.5 x 10 39 W 


2.5 x 10 39 W 


4.2 x 10 38 W 


^outcr 


2kpc 


2kpc 


3kpc 


lkpc 


lkpc 


lkpc 


B T 


2.65 x 10~ 6 T 


7.57 x HT 7 T 


1.02 x 10~ 6 T 


1.08 x 10~ 6 T 


5.93 x 10~ 7 T 


6.70 x 10~ 7 T 



Table 1. This table shows the values of the physical parameters used in the ballistic model fits shown in Figure[3] We also include the value of the magnetic 
field strength in the plasma rest frame at the transition region for the fits (Bt)- 



spectrum of J0531 if the transition region occurs at 10 R s as in the 
jet of M87. 

We have also considered an adiabatic model in which the jet 
is in equipartition after the transition region but with some electron 
acceleration occuring at smaller distances along the jet, in order 
to produce the gamma-rays and help reduce the bulk Lorentz fac- 
tor. In this case we find that accelerating electrons within the dusty 
torus tends to produce too much synchrotron emission due to the 
relatively high energy density in the magnetic field compared to 
the external photon field (both measured in the plasma rest frame) 
resulting in the overproduction of synchrotron emission. Acceler- 
ating electrons in the BLR also suffers from a similar problem un- 
less the location of the acceleration is fine-tuned to be close to the 
outer radius of the BLR. This is because the energy density in the 
BLR is relatively constant with distance x, whilst Ub oc x~ 2 . So 
the Compton-dominance, which is approximately proportional to 
Uj/Ub, increases towards the outer radius of the BLR. In an ac- 
celerating jet model the plasma has a lower bulk Lorentz factor at 
small distances within parabolic region and BLR, so the emission 



is not so strongly Doppler-boosted as that from further along the 
jet close to or beyond the transition region. We do not find that ac- 
celerating electrons within the BLR in the quiescent state is able 
to reproduce the observed data well and requires fine-tuning of the 
acceleration region within the BLR. For these reasons we do not 
favour this scenario in order to explain the quiescent spectra. 



6.3 Inverse-Compton scattering of external photons 

Our model uses a jet geometry based on radio observations of the 
jet in M87 scaled linearly with black hole mass. For large black 
hole masses (> 10 9 Mq) the distance of the transition region of 
the jet, where the jet is expected to be close to equipartition, is 
> 10 pc. This is consistent with our finding that BLR photons 
do not contribute significantly to the inverse-Compton emission of 
these Compton-dominant blazars since the BLR region is expected 
to be < 1 pc ( Kaspi et al. (2005 )). This means that the synchrotron 
bright transition region of our jet is outside the BLR and so the 
BLR photons will illuminate the jet from behind and be Doppler- 
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Figure 5. This figure shows the results of fitting the adiabatic jet model to the SEDs of J0349, J0457, J0531, J0730, J1504 and J1522. The model fits the 
observations well across all wavelengths, although for some of the blazars the model requires very high bulk Lorentz factors (7buik = 100 for J1522). We 
find that the inverse-Compton emission of all the blazars is well described using NLR seed photons with a contribution to the x-ray emission from scattering of 
CMB photons. 
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Parameter 


JUJ49.0 —2102 


J045 7.1— 262b 


J05il.U+iool 


Tf\T)f\ A 1 1 in 

J0730.4— 114z 


J 1504.4+1030 


J l52z.2+ol4o 


Wj 


8.0 X 10 W 


8.0 X 10 W 


1.4 X 10 W 


3.0 X 10 W 


2.0 X 10 W 


8.0 X 10 W 


L 


1 X 10 m 


2 x 1U m 


2 X 10 m 


1 X 10 m 


2 x 1.0 m 


2 X 10 m 


-^rnin 


5.1 1 Mev 


5.11 Mev 


10.2 Mev 


15.3 JvleV 


1U.2 JvleV 


15.3 JvleV 


^max 


1.28 GeV 


2.03 GeV 


1.14 GeV 


4.55 GeV 


2.61 GeV 


1.09 GeV 


a. 


1.3 


1.6 


1.55 


1.7 


1.9 


1.5 


nf 
O 

opening 


3° 


3° 


3° 


3° 


3° 


3° 


^observe 


1.25° 


1.0° 


0.8° 


1.0° 


0.8° 


0.6° 


70 


4 


4 


4 


4 


4 


4 


7max 


25 


50 


40 


50 


60 


100 


Tmin 


10 


30 


30 


30 


30 


30 


M 


8.0 x 10 9 M Q 


2.0 X 10 1U M Q 


4.2 x 1O 1U M 


2.0 x 1O 1U M 


1.5 x 10 1U M Q 


3.0 x 1O 1U M 


^acc 


5.0 x 10 4() W 


1.7 x 10 39 W 


6.8 x 10 39 W 


2.5 x 10 39 W 


8.3 x 10 39 W 


1.7 x 10 39 W 


^outcr 


Ikpc 


Ikpc 


Ikpc 


Ikpc 


Ikpc 


Ikpc 




4.02 x 10" 6 T 


8.05 x 10" 7 T 


6.34 x 10" 7 T 


1.56 x 10" 6 T 


1.41 x 10" e T 


2.68 x 10~ 7 T 



Table 2. This table shows the values of the physical parameters used in the adiabatic equipartition model fits shown in Figure [4] We also include the value of 
the magnetic field strength in the plasma rest frame at the transition region for the fits (Bt)- 



deboosted. We expect the transition region of the jet to dominate 
the high energy emission because the plasma is in equipartition 
and has the highest bulk Lorentz factor at this point giving the 
largest Doppler-boosting of the emission. From a reversal of the 
minimum energy argument (Burbidgej we know that for 

a given total energy, a plasma emits the most synchrotron power 
when the plasma is close to equipartition, which is found to be near 
to the transition region in jet simulations where the jet transitions 
from parabolic to conical ( McKinney ( 2006 1 and Komissarov et al. 
prjQ9l l). 

We show the different components of the inverse-Compton 
emission for the ballistic jet fit to blazar J 1504 in Figure [4] The 
figure shows that the inverse-Compton emission is dominated by 
scattering of CMB photons with a small contribution from scatter- 
ing NLR photons at high energies. For a jet with a transition region 
outside the BLR and dusty torus we expect CMB and NLR photons 
to dominate the external photon field (see Figure 4 in Paper II). At 
low redshifts we expect the energy density in starlight at the cen- 
tre of an elliptical galaxy to dominate CMB photons (Hardcastle 
|& Croston"|j201 1} ). However, these Compton-dominant objects are 
observed at relatively high redshifts so the CMB energy density is 
larger since pcmb oc (1 + z) 4 . So we expect the CMB to be the 
dominant scattered photon field several kiloparsecs from the central 
black hole. This is interesting since Compton-dominant blazars are 
observed predominantly at high redshift, which is exactly what we 
would expect if their inverse-Compton emission were due to scatter- 
ing CMB photons. Our model predicts that as we go to low redshifts 
where the energy density of the CMB decreases, high power, high 
bulk Lorentz factor blazars should become less Compton-dominant 
than an equivalent blazar at high redshift. We wish to investigate 
whether this prediction for the Compton-dominance of low redshift 
blazars holds in future work. 

In order for an electron distribution to produce a Compton- 
dominant spectrum from SSC emission we require that the syn- 
chrotron energy density is so high that a "Compton-catastrophe" 
occurs. This happens when electrons emit so much synchrotron ra- 
diation that the synchrotron radiation field is dense enough that scat- 
tering of photons is frequent and electrons lose more energy through 
inverse-Compton scattering than synchrotron emission. In a jet with 
a fixed total energy and equipartition fraction we require a small jet 
radius at the transition region in order for the synchrotron radiation 
field to be dense enough for a "Compton-catastrophe" to occur. We 
find that in such small radius/large magnetic fields the synchrotron 



radiation emitted by the highest energy electrons is much higher 
than the observed energy of the synchrotron peak frequency and so 
we are unable to fit to the observed spectra using SSC emission. 

6.4 Acceleration of electrons along the jet 

We find that the electron distribution power law indices are lower 
than expected from parallel shock acceleration (a — 2) for these 
Compton-dominant blazars. This suggests acceleration by oblique 
shocks pi ell et aT7| ( [20 1 1| | and Summerlin & Baring (20121), possi- 
bly a stationary recollimation shock at the transition region or accel- 
eration by magnetic reconnection (Zenitani & Hoshino (2001 1 and 
Zenitani & Hoshino (2007)). We see clear evidence for this in the 
spectrum of the blazar J053 1 where the synchrotron emission starts 
self-absorbed and close to flat and retains a steep gradient when 
the spectrum becomes optically thin (see Figure|5]5). This optically 
thin spectrum requires a steep electron spectral index to fit to the 
observed emission. In all these objects we also see that the inverse- 
Compton emission rises steeply from low energies to its peak. This 
steep peak is also consistent with a low electron spectral index. This 
is an interesting result since we found in Paper I that the spectrum of 
BL Lacertae was consistent with an electron spectral index a — 2, 
compatible with a parallel shock. We wish to investigate whether 
other BL Lac type objects require spectral indices that are higher 
than those of Compton-dominant objects. This could imply a dif- 
ference in the electron acceleration mechanism in different blazars 
possibly due to the different powers and bulk Lorentz factors of 
their jets. 

6.5 Compton drag 

In this investigation we have not explicitly calculated the effect of 
Compton drag on the jet plasma since we do not know the addi- 
tional amount of deceleration caused by the entrainment of matter 
surrounding the jet. We have used the simple prescription for the 
deceleration of the jet set out in Equation 5 of Paper II and fitted 
the deceleration parameters for a spectrum. We find for both the 
ballistic and adiabatic fits the majority of the jet's initial energy is 
not radiated away through anisotropic scattering of external pho- 
tons and so the effect of Compton drag on the jet plasma is com- 
parable to or smaller than half of the total amount of deceleration 
in the fits (A7b u ik < 0.27 max ). In the case of BL Lacs whose 
emission is dominated by synchrotron and SSC emission we do not 
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expect Compton drag to be a significant contribution to the jet de- 
celeration. This is because synchrotron and SSC emission from an 
isotropic electron distribution in a small scale tangled field are emit- 
ted isotropically in the plasma rest frame and so do not result in 
deceleration of the jet. This means that we are justified in not ex- 
plicitly calculating the deceleration due to Compton drag in these 
Compton-dominant FSRQs or BL Lacs since the contribution of 
Compton drag to the total deceleration of the jet is comparable to 
or smaller than that of entrainment. 



6.6 The blazar sequence 

We find that the slowly decelerating conical region, which is respon- 
sible for the radio emission, produces a nearly flat radio spectrum 
which is consistent with observations. The frequency at which the 
observed synchrotron spectrum transitions from optically thick to 
optically thin emission is governed chiefly by the black hole mass 
which sets the radius of the bright transition region (the base of the 
conical region) and largely determines the magnetic field strength at 
this point. The transition from optically thick to optically thin syn- 
chrotron emission at these comparatively low energies (as shown 
in the spectra of J0457 and J0531) implies large black hole masses 
for these blazars. This is because a transition to optically thin syn- 
chrotron emission at low frequencies implies a low magnetic field 
strength at the base of the conical region and therefore a large radius 
as shown in Equation[8] 

The base of the conical region in our model is in equiparti- 
tion and so is responsible for producing the highest energy optically 
thin synchrotron emission. The peak frequency of the synchrotron 
emission is then governed by the Doppler factor, the magnetic field 
strength at the base of the conical region and the maximum elec- 
tron energy. We find that the low synchrotron peak frequency of 
these Compton-dominant blazars is due to their large black hole 
masses which result in the transition region occuring at large dis- 
tances from the black hole and with large radii. These large radii 
result in relatively low magnetic field strengths at the bright tran- 
sition region which produce the low peak frequency synchrotron 
emission. We therefore expect that low powered BL Lac objects 
with high synchrotron peak frequencies should have lower black 
hole masses than low synchrotron peak frequency objects such as 
Compton-dominant objects. 

This can be illustrated through the following simple calcu- 
lation. We assume that the base of the conical region at a dis- 
tance 10 5 r s from the central black hole with a radius 2000r s is in 
equipartition and has the maximum bulk Lorentz factor 7 max , so we 
expect it to dominate the optically thin synchrotron emission. For 
an electron distribution with a maximum Lorentz factor 7 e and an 
initial kinetic Luminosity Wj using the formula for the synchrotron 
critical frequency (Equations [3] and |4j and Equation[l3]we find 



B 2/i 3 ~ 2c ' 



C M 

Rt = 2000^^, (13) 



^pcak 



: ppC/'y" max 



3<5d. 



47r(l + 2)m|c 4 y 47, 



3W>p 
c(2000GM/c 



fpcak « W] /2 /M_ 



(14) 



(15) 



where a subscript T denotes the value of a quantity at the tran- 
sition region, Mbh is the central black hole mass, j3 is the velocity 
of the jet divided by c, 9 h s is the observation angle to the jet axis, 



z is the redshift and <5do PP is the doppler factor for the jet. We have 
assumed that half the initial jet power is contained in magnetic field 
at the transition region where the plasma is in equipartition. We 
have also assumed the plasma is described by a relativistic perfect 
fluid in its rest frame. From this we see that for objects which differ 
only by black hole mass z^ pC ak oc 1/Mbh- If we further assume 
that in general the jet power is a fixed fraction of the Eddington lu- 

— 1/2 

minosity then Wj oc Mbh and so t'pcak oc M S£r . If blazar jets 
are reasonably well described by a scaled M87 geometry this could 
even be used as a rough estimate of the black hole mass of a blazar 
based on its synchrotron peak frequency. 

It would be an exciting result if the black hole mass of a blazar 
largely determines the peak frequency of its synchrotron emission. 
It is perhaps not surprising if in general FSRQs have higher black 
hole masses than BL Lacs. If jet power is fundamentally related 
to accretion, as currently believed, then we would expect powerful 
FSRQs to be accreting at a higher rate, and therefore to have larger 
black hole masses than lower power BL Lacs. In this scenario we 
would expect to find a trend like the blazar sequence where power- 
ful objects are Compton-dominant and have lower peak frequencies 
of emission than less powerful BL Lac type objects with smaller 
jet radii (due to smaller black hole masses) and correspondingly 
higher magnetic field strengths. We interpret this as due to the bulk 
Lorentz factor of the jet, redshift and the black hole mass (which de- 
termines the size of the jet). If we assume a jet geometry similar to 
that observed in M87 we expect the jet to transition from parabolic 
to conical at 10° R a . The jet radius at this point is a function only 
of black hole mass (if we assume the M87 geometry scaled linearly 
with black hole mass) and so the magnetic field strength and syn- 
chrotron peak frequency decrease with increasing black hole mass. 
For BL Lac type objects with low black hole masses the radius of 
the jet at the point where it transitions to conical and is in equipar- 
tition is smaller than in a more powerful source with a higher black 
hole mass and this results in the increased magnetic field strength. 
An empirical correlation between jet power and black hole mass has 
been found by several authors by comparing observed radio lumi- 
nosity against black hole mass for AGN JLacy et al.|pOOT] , |McLurel 
|& Jarv is (20041 and Liu et al. (2006}), however, this correlation has 
been disputed ( |Oshlack et al] p002l >, [Woo & Urry| {20021 and [Ho] 
p002l ». 

If we assume that all AGN jets have approximately the same 
geometry as the M87 jet we can derive a simple relationship for 
the peak synchrotron frequency in blazars. In fitting to the optically 
thick to thin synchrotron break of these blazars we are essentially 
fitting the size of the jet radius where the synchrotron emission is 
brightest. In our model we then associate this emitting region with 
the region where the jet transitions from parabolic to conical. The 
determination of this radius from fitting to a spectrum is a rigorous 
result which holds independently of the particular M87 jet geome- 
try we have chosen. In this work our assumption that the geometry 
of the jet scales linearly with black hole mass is based on the expec- 
tation that the radius of the jet close to the black hole scales with the 
size of the black hole event horizon which depends linearly on the 
black hole mass. Of course it is likely that the geometry of individ- 
ual jets depends on other quantities such as jet power, bulk Lorentz 
factor and the environment. We do not yet know from observations 
or simulations exactly how these factors affect the jet geometry so 
we assume a simple scaling with black hole mass. It is important 
to note that if the inner parabolic geometry of jets is not well de- 
scribed by this simple model this will mainly affect the inferred 
black hole mass of our fits and not the other physical parameters 
of the model. This is the first investigation which has attempted to 
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constrain the physical parameters of a sample of blazars using a 
realistic, extended model. 

The trend of decreasing peak frequency of inverse-Compton 
emission with jet power can also be explained using this model. 
Low powered BL Lac type objects with low black hole masses have 
relatively high magnetic field strengths at the base of the conical 
section. This means that the synchrotron emission will overwhelm 
CMB photons (especially since BL Lacs are observed at lower red- 
shifts than FSRQs and have lower bulk Lorentz factors) and the 
inverse-Compton emission will be due primarily to scattering of 
relatively high frequency synchrotron photons. This results in an 
inverse-Compton peak occurring at a higher frequency in BL Lacs 
than in Compton-dominant objects. This is because the inverse- 
Compton peak frequency is proportional to the peak frequency of 
the seed photons ^ pcak oc 7 2 ^ sco d (for i? max » /i^ pca k) and syn- 
chrotron photons have a higher peak frequency than CMB photons. 

For increasing black hole mass, the magnetic field strength de- 
creases at the base of the conical section. The synchrotron peak fre- 
quency decreases and this means the SSC peak frequency will also 
decrease. For higher powered, high bulk Lorentz factor jets with 
lower magnetic field strengths Doppler-boosted CMB photons will 
now dominate the photon distribution within the jet at large dis- 
tances and electrons will lose most of their energy scattering CMB 
photons instead of emitting synchrotron photons. This produces the 
Compton-dominance of these objects and the low inverse-Compton 
peak frequency is due to the comparatively low CMB frequency. 

In this paper we have not attempted to explain the short 
timescale variability in Compton-dominant blazars but instead the 
quiescent emission. The models we have investigated, which fit well 
to the quiescent emission, have electron acceleration occuring at 
large distances along the jet where the light-crossing time of the jet 
and electron cooling times are both much larger than the short flar- 
ing timescales observed (hours/minutes Aharo man et al.| |2007[ >). 
These flares could, however, be produced by acceleration of elec- 
trons in the parabolic region of the jet where the radius of the jet 
is smaller, the magnetic field strength is higher and the jet is still 
within the BLR/dusty torus. This means that short timescale vari- 
ability is compatible with our jet model. 

In this investigation we have shown that a ballistic jet model 
which becomes particle dominated beyond a transition region lo- 
cated > lOpc fits the multi wavelength data of the six most 
Compton-dominant blazars (after PKS0227) from |Abdo et al.| 
( |2010[ >. We have demonstrated that using a realistic extended jet 
model allows us to gain constraints on the jet by reproducing radio 
observations. In the next paper in the series we wish to investigate 
if our model is able to reproduce the spectra of BL Lac type blazars 
and whether our idea of a blazar sequence based on bulk Lorentz 
factor and black hole mass holds for a sample of BL Lacs. 



7 CONCLUSION 

In this paper we have used the realistic extended jet model from 
Paper II to fit to the simultaneous multi-wavelength spectra of 
six Compton-dominant blazars. Our model includes a magnetically 
dominated accelerating parabolic base transitioning to a slowly de- 
celerating conical jet with a jet geometry set by radio observations 
of M87 and consistent with simulations and theory. Comparing dif- 
ferent models of the conical section to the radio data, we find that 
the radio slope is well matched by both an adiabatic equipartition 
jet and a ballistic jet which starts in equipartition, at the transition 
from parabolic to conical, and becomes particle dominated further 



along the jet. We use the optically thick to thin synchrotron break to 
constrain the radius of the jet where the jet first comes into equipar- 
tition. We find that this requires a transition region located outside 
the dusty torus at large distances along the jet > lOpc, consistent 
with a transition region occuring at 10° R s as in the jet of M87. We 
calculate an analytic expression for the radius at which a plasma 
becomes optically thick to an observed frequency and use this to 
confirm our numerical results. We then use this formula to calculate 
the expected frequency core-shift relations for the jet models under 
consideration. We find that for the adiabatic equipartition conical 
jet the core-shift relation, x oc v~ k , has a value k — 0.69, lower 
than suggested by observations. 

We find that the ballistic particle dominated jet fits very well 
to the observations across all wavelengths including radio observa- 
tions whilst the adiabatic equipartition jet requires very high bulk 
Lorentz factors to reproduce the Compton-dominance of some of 
the blazars 7b u ik > 50. This is the first time a realistic, extended 
jet model has been used to fit to a sample of blazars. We find that 
the fits to all these blazars require high power (> 10 39 W), high 
bulk Lorentz factor (> 20) jets observed close to the line of sight 
(< 2°) consistent with the blazar sequence and unification and in 
agreement with the results from Paper II. We find that the inverse- 
Compton emission of these objects is well fitted by scattering of 
CMB seed photons due to the redshift dependence of the CMB and 
strong Doppler-boosting. We find that all these blazars have large 
inferred black hole masses (> 1O 9 M0) using a jet geometry set 
by radio observations of M87 scaled linearly with black hole mass. 
For the blazar J0531 we find that the black hole mass inferred from 
fitting the jet model to the optically thick to thin synchrotron break 
agrees with the black hole mass estimated independently from the 
prominent accretion disc spectrum only if the transition region oc- 
curs at 10 J -R S as in M87. For these large black hole masses, con- 
sistent with massive elliptical galaxies, we find that the brightest 
part of the jet is located at 10°-R S , outside the broad line region and 
dusty torus. 

We find a relatively low value for the electron distribution 
power law index from fitting to the steep optically thin synchrotron 
emission and inverse-Compton peak (a < 1.9), implying accel- 
eration by oblique shocks or magnetic reconnection. We postulate 
a new interpretation of the blazar sequence in terms of increasing 
black hole mass, jet power and bulk Lorentz factor. Higher power 
blazars have larger black hole masses and larger radius jets with 
lower magnetic fields than low power BL Lacs. This results in 
high power blazars having lower synchrotron peak frequencies and 
less SSC emission than low power blazars with higher magnetic 
fields. High power blazars also have large bulk Lorentz factors and 
so scatter high-redshift CMB photons at large distances resulting 
in the low inverse-Compton peak frequency of Compton-dominant 
blazars. 

These results illustrate the importance of using a realistic, 
physically motivated model for blazar jets to calculate their spec- 
tra and estimate their physical properties. In a future paper we will 
test these concepts further by fitting to the spectra of a sample of 
BL Lac type objects. 
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